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Abstract. 

We show that the quantum fideUty is accessible to cold atom experiments for a 
large class of evolutions in periodical potentials, properly taking into account the 
experimental initial conditions of the atomic ensemble. We prove analytically that, 
1/^ , at the fundamental quantum resonances of the Atom-Optics Kicked Rotor, the fidelity 

saturates at a constant, time-independent value after a small number of kicks. The 
. latter saturation arises from the bulk of the atomic ensemble, whilst for the resonantly 

I accelerated atoms the fidelity is predicted to decay slowly according to a power law. 

o 



PACS numbers: 42.50.Vk,03.75.Be,32.80.Qk,05.60.Gg 
1. Introduction 



^1 One of the most intriguing quantum phenomena is the coherent evolution and control 

of a superposition of state vectors, which is at the heart of quantum parallelism in 
^ ' quantum computation schemes. One way to directly probe coherent quantum evolution 
' is offered by the fidelity, which correlates the time evolution of a reference system with 

2 

its slightly perturbed dynamics through F{t) = {■ipQ\e^-^'^''^^^e~^^^^^\-ipo) ■ Here, H 

and = H + AV represent the reference Hamiltonian and its perturbed variant, 
respectively, which generate the evolution, starting from the same initial state {ipo). F{t) 
was introduced as a measure for the stability of quantum motion with respect to changes 
in some control parameter A jTj. The fidelity has recently evoked considerable interest as 
an alternative route for the study of the effects of perturbations on the coherent evolution 
of quantum systems, in particular in the context of quantum-classical correspondence j2] 
and of quantum information processing jH] . As for the former, the finite size of h makes 
the classical definition of chaos through the exponential divergence of initially nearby 
trajectories meaningless on the quantum level. However, the sensitivity to changes in 
A prevails both in the classical as well as in the quantum realm, and is dynamically 
quantified by F(t). 

Whilst previous studies focused on the limit of small effective Planck constant 
[3 E] , we present in this Letter an analytic result which is valid (a) in the deep quantum 
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regime, and (b) for arbitrary values of the control parameter A. We study the fidelity 
for cold, periodically kicked atoms, which mimic a paradigmatic model of quantum 
chaos theory, the (5-Kicked Rotor (KR) In this context, we show that bona fide 
fidelity is readily accessible, yet has never been measured in an experiment such as 
reported in Specifically, we derive an analytic expression for the fidelity at the 
fundamental quantum resonances of the KR. At these resonances, a subclass of the 
experimental atomic ensemble [3 |Hl El experiences an enhanced acceleration by the 
kicking field, at precisely defined kicking periods [H]. We show that this resonantly 
driven subclass induces a slow algebraic decay of fidelity, whilst the bulk of the initial 
atomic ensemble produces the novel phenomenon of a perfect saturation of fidelity at a 
finite and experimentally easily detectable level. 

Quantum resonances are of purely quantum origin, although it has been recently 
shown that they can be mapped onto a pseudo-classical problem. In this sense, we 
study the quantum fidelity in the deep quantum regime, far from the semiclassical limit 
of large actions, but close to the pseudo- classical limit. The latter anchors the near- 
resonant quantum dynamics to elliptic, stable resonance islands in the phase space of a 
suitably defined classical map, precisely alike the usual semiclassical limit. However, the 
pseudo classical limit is reached at fixed and finite h [10^. Notwithstanding this formal 
complication, the correspondence between (pseudo) classically stable dynamics and near- 
resonant quantum evolution immediately suggests the remarkable robustness of the 
quantum dynamics at resonance, originally believed to be very sensitive to perturbations 
(see fl] for a further discussion). Hence, quantum resonances provide a novel and robust 
tool of quantum control, and fall within a larger class of particularly robust modes 
originating from the dynamics of strongly perturbed quantum systems [TIH I12j. This 
robustness manifests indeed in the temporal evolution of fidelity as shown in the sequel. 

2. The Atom-Optics Kicked Rotor and Quantum Resonant Motion 

The experimental realization [71 |H] of the KR uses cold atoms periodically kicked by 
an optical lattice. The Hamiltonian that generates the evolution of the atomic wave 
function in its centre-of-mass degree of freedom is given in the following dimensionless 
form [T^ 



where P is the momentum operator, and X the atom's position along the pulsed optical 
lattice. The two parameters that control the dynamics are the kicking period r, which 
takes the role of an effective Planck constant in our units 0, and the kick strength 
k. The Hamiltonian describes the Atom-Optics Kicked Rotor (AOKR) with the 
atoms moving along a line rather than on a circle (as in the standard KR model). 
(H} induces the characteristic properties of the quantum KR, which sensitively depend 
on the parameter r CD]. If t = Anr/q, with r, q integers, then the kicking period is 




(1) 
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rationally related to the propagation time of the kicked atoms across the lattice constant 
(which is 27r in our units), and quantum transport is enhanced. Otherwise classically 
diffusive energy growth is suppressed by dynamical localisation In the following, 
we derive analytical results for the fidelity, for the experimentally relevant quantum 
resonances at r = 27ii e M) P], i.e., for large effective Planck constants r. 

An atom subject to a potential periodic in space and time is described by a 
wave packet composed of 27r-periodic Bloch states, ip/3{x + 27r) = ip/3{x), with the 
quasimomentum /3. In our units, the latter is given by the fractional part of the 
momentum p = n + /?, n G N. /5 is a constant of motion, so the different Bloch states 
evolve independently of one another, and their momenta only change by integers. At 
T = 27r£, a special behaviour occurs only for a specific, small subclass of quasimomenta f3. 
In much the same way as spatial periodicity enforces ballistic motion in physical space, 
the one-period evolution operator is periodic in momentum space, enforcing ballistic 
propagation of the corresponding /3-states in momentum space. As shown in ^U], during 
one period r, \ipi3) evolves into Up \1p13), with 

Up = e-^^™«W e-^5(^+^)', (2) 

where 9 = Xmod(27r), and J\f = — i d/d9 is the angular momentum operator. For 
instance, at r = An, and (3 = 0, 1/2, or at r = 27r and (3 = 1/2, the second factor of 
Up equals the identity (apart from a global phase for /3 7^ 0), and the energy of the 
corresponding /3-state grows ballistically, i.e., quadratically with the number N of kicks 
0. For general values of /5, but r = 27ii, i integer, the state after the iVth kick is, in 
the momentum representation [TO] 

/•2iT rn 

{n \U^tpp) = e'"^'^^^^^) / e-^^^-'^^I^^I^^^^V/s {9 - m - arg(W^jv)) , (3) 

Jo V 27r 

where = W]\f{^) := X^fL^^ exp(— is^). If the initial state of the atom is a plane wave 
of momentum po = tiq+P, then ^ takes the constant value ^0 = vr£(2/3— 1). If ^0 7^ 0, then 
\Wn\ = |sin(A^^o/2)/sin(^o/2)|, else if ^0 = then \Wn\ = N, and the average kinetic 
energy increases like The latter occurs at the resonant values f3 = 1/2 -\-n/i 

mod(l), n = 0, !,..,£— 1, which support ballistic motion ^U]- The remaining Bloch 
components of the original wave packet, with f3 not belonging to the resonant subclass, 
undergo a quasiperiodic energy exchange with the kicking field, leading to a finite spread 
of the associated momentum distribution for all times ^U]. How does a perturbation 
in the kick strength k affect such quantum dynamics at quantum resonance - brought 
about by the abovementioned complete phase revival in-between two successive kicks? 
The answer to this question will be given by computing the fidelity of the resonant time 
evolution, using eq. for two different values of k. Before doing so, we show how to 
extract the fidelity from observables directly measurable in the experiment 0. 
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Figure 1. Experimental setup of ,6 . Each atom is prepared in the upper state (left), 
then a 7r/2 Ramsey pulse is applied. The dynamics of the coherent superposition of 
the upper (where k = ^2) and lower {k = ki) level lasts for N kicks, before a second 
it/ 2 Ramsey pulse is applied (with time delay T^a with respect to the first one). In 
principal, any type of dynamics may occur in-between the two Ramsey pulses, but 
our analytical results restrict here to the AOKR, as an example for a coherent atomic 
evolution in a periodical optical lattice. 



3. Experimental Measurement of Fidelity 

The experiment uses two sublevels of the atomic ground state of the electronic degree 
of freedom, which couple differently to the kick potential. Before the kicks are applied, 
the sublevels are equally populated, and then evolve coherently. The splitting of the 
population is performed by a Ramsey- type interferometer technique jHl the 
complete time evolution of the atoms is sketched in Fig. ^ The internal states are 
represented by V^^^"*), where the upper index denotes the corresponding electronic 
level of the atom. Ei,E2 are the energies of the two levels. The Ramsey pulses are 
tuned to the electronic transition between these two levels, i.e., they have a frequency 
cura — {E2 — El) /h. For simplicity, we assume instantaneous Ramsey pulses which occur 
immediately before the first and after the last kick, respectively. These pulses couple the 
two states whilst they evolve coherently, but independently, under the kicking potential. 
Let W^^'* and W^^'* be the propagators (j2I) in the upper and in the lower atomic level, 
respectively, and (f) the relative phase of the two Ramsey pulses. For an initial state where 
all the population resides in one of the electronic sublevels, i.e., V'^^'' = ^ip^^^ = ipis, 
the states after the final Ramsey pulse are [TT] 

4\N,S = V-.- ((M»')"e-*" + (Wf )") V, 

with the phase (/)ra = {E2 — Ei)N — cj). Projection onto a momentum eigenstate with 
momentum n gives the momentum distribution, and we obtain the total probability in 
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N 



(4) 



state ipp'^ upon summation over all n: 

The observable given by eq. can be accessed experimentally by measuring the atomic 
momentum distribution. The measurement, however, provides only an average over all 
(5 in the initial ensemble [Oj. Thus, when the initial pure state ipp (of the atomic centre- 
of-mass motion) is replaced by a mixture which is described by the statistical density 



operator p = d(3 f{(3)\i)p){i)p\ 

1 
2 

We defined [II] 



the final, measurable population in state V'L^^ is 



Pl(iV,0aA) = ^ {l + COS(0")} 



COS 



d(3f{P)Re 







\/F{N) 

where, at fixed time A^, 0" differs from 
F{N) 



'/3 



N 



N 



(5) 



(6) 



N 



as well as from by a constant shift, and 

2 



/(2) 



N 



(7) 



corresponds to the fidelity for an incoherent ensemble of /^-states. From this analysis, 
we see that the fidelity is equal to the squared difference between the maximum and the 
minimum values taken by Pi{N, 0), while the original phase varies in [0, 2%) pT| ITK]. 

Hence, the fidelity is accessible as the visibility of the integrated momentum 
distribution as a function of the Ramsey phase by the experimental setup of 
Fig. ^ One must, however, explicitly take into account the scan of the Ramsey phase, 
in contrast to what has been dubbed "fidelity" in [Hj. Moreover, as long as the evolution 
operators are known and the atoms evolve in a spatially periodical potential, eqs. (^HZj) 
remain valid and are not restricted to the here analysed AOKR. 



4. Analytical Results and Discussions 

For the AOKR experiments [71 |H] , with an initially Gaussian momentum distribution 
with root mean square > 1, the distribution of quasimomenta /o(/3) is nearly uniform 
|lUj . Since /? is a constant of motion, the fidelity is 



FiN) 



p +00 



To compute F{N) we must calculate the overlap of two wave functions Q for two values 
of k, and sum over all integer momenta n. For the case when the state of the atom is a 
plane wave with momentum Po = Uq + /?, we have |'0/3(^ ~ ~ arg(1^7v))P = l/(27r), 
leading to the result for the fidelity for one /3-state: Fj3{N) = Jq (|VFiv|(A;i — k2)). For 
Ak = ki — k2 = 0, Fi3 = Fi3{N) = 1. For resonant /3 we have \ Wn\ = N, and we can use 
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the asymptotic expansion formula [9.2.1] from ^B] for the 0th order Bessel function Jq. 
The fidehty for large times at quantum resonance is then 

Fa (N) ~ cos^ ( AkN - -) . (8) 

This shows that the fidelity falls off as a power law oc 1/A^ for the resonant /^-states. 
Another particular case is given for (3 = 0: Fg=o(-^) = ./q (| sin(A^7r/2)|AA;), which at 
= 1 depends just on Ak. Then the fidelity can be made arbitrarily small when 
Ak ~ zq, where Zq ~ 2.40 is the first zero of Jq. In turn, the fidelity can be maximised 
by maximising Jq, what provides an option for coherent control when quasimomentum 
can be precisely controlled in experiments with Bose-Einstein condensates 0. 

To compute F{N) for fo{/3) = 1, what approximates very well the experimental 
situation in jSJ^ni! we perform the average over the different quasimomenta, i.e. 

r dx 

dpJo{\WN{0\^k) = / —Jo{Aksm{Nx)csc{x)) 



271 



N-l 



f^^j7Jo{^ksin{x)csc{xN-^ + 2TxrN-^)) . (9) 

-'^ ■r=0 

We used the definition of Wn^C) as given after eq. ©, and the periodicity of the sine in 
the argument of the Bessel function. In the limit when N oo and 2'Kr/N — > a, the 
sum over r approximates the integral over a, and (jHl) converges to 



d(3Jo{\WN{C)\^k) 7K-^ I dx I da Jq {Ak sm{x) csc{a)) 

(27r)^ 



With the help of [11.4.7] from [16^, this leads us to the final result for the asymptotic 
value F* of the fidelity HH 



1 / Z"^" . fAkcsciaV ' ^ 



^•<^*)"(wU ■ (1°) 

The dependence on A A; of the formula ()1()|1 is shown in Fig. |21 together with numerical 
computations. F*{Ak) oscillates quasiperiodically rather than to drop monotonously 
with Ak. Note that the fidelity always saturates at exact resonance r = 27r£ (£ G N), and 
the saturation time and level depend on Ak, as illustrated in the inset. We interprete our 
analytical result (fTO|l by referring to the wave function in eq. (jH)) : the perturbation enters 
in the expression for the fidelity in the form of phase factors exp (— i A/c|VI/Ar| cos(^)). 
For non-resonant (3, \Wn\ — > in the limit N — > oo, and hence these phases eventually 
tend to unity, i.e., the perturbation completely looses its effect on the relative evolution 
of the two wave packets under comparison. Note that a fidelity plateau, somewhat 
reminiscent of the saturation observed here, was recently predicted for some model 
systems j3]. However, this strictly semiclassical result was obtained in the perturbative 
limit, and the plateau observed in ji] only represents a transient effect, followed by 
the subsequent decay of fidelity. By contrast, our formula pO|l is a result valid for all 
(sufficiently large) N. 
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Figure 2. Asymptotical formula for the fidelity for iV — > cx) from eq. H10|l (solid 
line), compared with numerical data (squares) obtained by evolving ensembles of 10^ 
/3— rotors with an initially uniform momentum distribution in [0,1), at r = 27r and 
N — 50. Inset: numerically computed fidelity vs. NAk for r = 2n, ki — O.Stt, and 
fixed Afe = 0.6283 (crosses), 1.885 (diamonds). The position of the minima corresponds 
to the time S.SSr/Afc. The fidelity saturates for times > 20 at Afc-dependent, 
constant values, indicated by the horizontal lines. Data for finite detunings e = 0.025 
(dashed) and 0.1 (dot-dashed) is shown for Afc — 0.6283. 



The asymptotic saturation value F* depends on the difference Ak in the kick 
strength and, in addition, on the distribution of the initial conditions /o(/3): the resonant 
f5 - which behave according to (jH)) - determine the initial decay and the consecutive rise 
of the fidehty F{N) shown in the inset of Fig. |21 Their contribution quickly decreases, 
since the speed of the wings of the two wave packets to be compared is proportional to 
the two different k jHlEI]. After this initial stage, the bulk of non-resonant (3 takes over 
and dominates then the fidehty, leading to rapid saturation. These competing effects 
give rise to the damped oscillations around the asymptotic value of fidelity. At fixed Ak, 
F* is the lower, the more weight is given to quasimomenta which are close to resonant 
values. The initial decay of the fidelity is due to nearly resonant (3 values: the minima in 
the inset in Fig. |21 are found by differentiating the time- dependent fidelity with respect 
to Ak. Then, the first zero of the 1st order Bessel function Ji{z) (z = ~ 3.83) 
is responsible for the observed minima. Figure |21 also shows that strict saturation 
is destroyed by arbitrarily small detunings e = 27r£ — r from resonance (since then 
the quasiperiodic motion for the non- resonant (3 is destroyed). However, the temporal 
asymptotic decay of the fidelity depends continuously on e, and is very slow for small 
e, what leaves the here predicted saturation an experimentally robust observable (since 
the kicking period is the best controlled parameter in the experiment lEllZllHl)- 
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5. Conclusions 

In summary, we analytically derived the asymptotic behaviour of fidelity for kicked cold 
atoms at quantum resonance, far from the classical limit and for arbitrary perturbations 
in the kicking strength k. The predicted saturation sets in after a small number of 
kicks, at an experimentally observable level. Note, in particular, that F* > 50 % for 
Ak ~ 0{k), what clearly stresses the robustness of quantum resonant motion. The 
type of dynamics occurring in-between the two interferometric pulses (see Fig. is not 
restricted to the AOKR. Our derivation [see eqs. (jHITj)] applies for arbitrary evolutions 
of cold atoms in a periodical potential, since it properly takes into account the different 
quasimomentum classes (corresponding to all the initial conditions) of the experimental 
ensemble. 
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